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I ■ concerns the asymptotic behaviour of the componentwise sample minima related to an elliptical distributions. Further, 
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Introduction 



Let R and S be two independent positive random variables. In this paper we consider the random scaling model 

W ^ RS, (1.1) 



, with W the scaled version of R and S e (0, 1) almost surely (= stands for equality of the distribution functions). In 
^ ' order to derive distributional properties of W we need to specify the distribution function of S; a tractable instance 
, with various applications is the tractable case that S" is a beta distributed random variable. 

00 

\0 ' 

In a financial or insurance framework, the random scaling model (jl.ip appears naturally with W the deflated risk 
■ arising from some loss or investment R which is independent from the random scaling/deflating factor S. Other 
' prominent applications in the literature concern modeling of network data (see e.g., D'Auria and Resnick (2006, 2008)); 
. random difference equations (see e.g., Mikosch and Konstantinides (2004), Denisov and Zwart (2007)); insurance and 
; finance applications (see e.g.. Tang and Tsitsiashvili (2003, 2004), Tang (2006, 2006), Piterbarg et al. (2009), Liu 
I> . and Tang (2010), Tang and Vernic (2010), Zhang (2010)); approximation of multivariate distributions (see e.g., 
I Hashorva (2007), Charpentier and Segers (2009), McNeil and Neslehova (2009), Balakrishnan and Hashorva (2010)). 
, A monograph treatment rich in applications and references is Galambos and Simonelli (2004). 
' 

■ - - ' When R is not directly observable, but the distribution function (df) of S is known, and W is observable, a natural 
question arising from (jl.ip is the recovery of the distribution function of R, or its estimation. Such a question arises 
for instance while estimating the true claim cost of a glass insurance coverage. Indeed, if i?^, i > 1 models the losses 
payed to claims reported from some glass coverage of a particular motor portfolio, the insurer is interested in the 
estimation of the true claim cost Wi . However, these costs are typically deflations of Ri , where the deflator Si explains 
the presence of fraud or other effects; in this setup Ri is not directly observable. 

In certain cases the df of the scaling random variable is known, or it can be estimated, which prompts the insurer to 
attempt to recover the df of the true losses. This is possible when the df of the random variable 14^ is a beta-product 
convolution, i.e., the scaling random variable S is beta distributed with positive parameters a, /3, see (jS.lip below. An 
interesting fact connected with beta-product convolutions is the characterisation of fc-monotonc functions, see Fakes 
and Navarro (2007), Balabdaoui and Wellner (2010). 

The principal aim of this paper is the investigation of the scaling model (jl.ip from both distributional and asymptotical 
interest extending some previous findings of Hashorva and Fakes (2010) and Hashorva et al. (2010). In the framework 
of beta random scaling we consider the inverse problems of derivation of the distribution (or the density) function 
of R when that of W is known, which can be solved by resorting to properties of the Weyl fractional-order integral 
operator. In the second part of the paper we deal with the lower tail asymptotic behaviour of W and R, dropping 
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specific distributional assumptions on S. In particular, we investigate the min-domain of attraction of the distribution 
function of W if that of R belongs to the min-domain of attraction of some univariate extreme value distribution 
function. For the beta-product convolution model we are able to derive further some converse asymptotic results. 

Distributional properties and results such as lower (upper) tail asymptotics of beta-product convolutions are of certain 
importance for insurance application when dealing for instance with the modeling of small (large) claims which are 
typically affected by some random inflation (deflation) factor. In fact, from the financial point of view, insurance 
companies do not suffer from small claims but from the large ones. However, understanding small claims is important 
for at least two reasons: a) claim handling is expansive even for zero-losses or very small ones, b) the choice of 
deductibles and the calculation of pure premiums can be significantly improved if the effect of inflation/deflation on 
small claims is adequately modeled. In finance, modeling of the effect of a defiator, which can practically ruin an 
investment, is very important. 

We present in this paper two applications: first we investigate the asymptotic behaviour of the componentwise minima 
of the absolute value of elliptical random vectors where we show that the it is attracted by some multivariate df with 
independent components, provided that the associated random radius a regularly varying (at 0) df. In the second 
application we consider the aggregation of two risks with polar representation (similar to that of elliptically distributed 
risks). Aggregation of risks is an important topic for insurance and finance, see the recent contributions Dhaene et 
al. (2008), Asmussen and Rojas-Nandaypa (2008), Embrechts and Puccetti (2008), Albrecher and Kortschak (2009), 
Geluk and Tang (2009), Mitra and Resnick (2009), Valdez et al. (2009), and Dengen et al. (2010). 

Outline of the rest of the paper: Preliminary results will be followed by Section 3 where we discuss focus on the main 
distributional properties underlying the beta random scaling model. 

Lower tail asymptotics for W and R related by is investigated in Section 4. Two applications in Section 5 

proceed the las section which contains some related results and the our proofs. 



2 Preliminaries 

In the sequel a, /? are two positive constants, and Ba^p denotes a beta random variable with density function 

r(a)r(/3) ^ ' ' ''^^^''>' 

where r(-) is the Euler gamma function. The distribution function of R will be denoted by H (abbreviate this as 
R ~ H) Ha,i3 is the df of W with stochastic representation with upper endpoint of ujuj G (0,oo]. It will be 

assumed that the lower endpoint of 77 is O (so H{jS) = 0). For our scaling model (jl.ip the df of W is said to be a 
product convolution distribution defined in terms of H and the df of 5*. When S is beta distributed with parameters 
a, 13 the relation between H and -ffa,/3, with H^^p the df of S is quite tractable due to the role of the Weyl fractional- 
order integral operator. We refer to J?a,/3 alternatively as a beta-product convolution. 

Next, we introduce the aforementioned operator acting on real- valued measurable functions h defined on (0,oo). For 
a given constant /3 € (0,oo) the Weyl fractional-order integral operator Ip is defined by 

1 1"°° 

mh){x) ^ Y\p)J iy-^f~'hiv)dy^ ^>o. 

Now, if for any £ > we have 

/•oo 

a;^~"'"|/i(a:)| dx < oo, 



which is abbreviated by G I^, then {Iph){x) is almost surely finite for all x € (0,oo). 
It follows easily that 

Ho.A^) = ^%^a;"(/^p_„-/3i?)(x), xeiO,Lj), with p,(x) = s e (0, oo) (2.1) 

r(a) 

showing the importance of the Weyl fractional-order integral operator in the setup of beta random scaling. 
When a = 1, then P{Bi,;3 > s} = (1 - s)'^, s e (0, 1). Hence simplifies to 



HcA^)= {l-x/yfdH{y), H^,p = l-H^,^ 

J X 
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which leads us to the introduction of the Weyl-Stiehjes fractional-order integral operator Jp^g with g : (0, oo) -^M a 
measurable weight function defined by 

1 f°° 

{Jp,gH){x) = — - / {y-xf-^g(jj)dH{y), xe(0,c.). 



With this notation we have 

Hi,p{x) = r(/3 + l)(J'^+i,p_,i/)(x), xe(0,c^). (2.2) 

When /3 = d ElN, then T{P){Jp,p^_^H) is the Williamson d-transform of H . which plays a crucial role in the analysis 
of Li-norm Dirichlet distributions (Fang et al. (1990)), and Archimedean copula (McNeil and Neslehova (2009)). 

For the derivation of the lower tail asymptotics of W we impose an assumption on R motivated by univariate extreme 
value theory. Specifically, we assume that R ^ H is regularly varying at with some index 7 G (0, 00), i.e., 

lim^^^ = Vie (0,00). (2.3) 

xio H[x) 

Alternatively, we write H £ RVy or R £ RV^ . Eq. (|2.3p is equivalent with the assertion is regularly varying at 
infinity with index —7, or the df of 1/R is in the max-domain of attraction of the Frechet df '^7(2;) = exp(— x"'''), a; > 0. 
See Resnick (1987), Bingham et al. (1987), Embrechts et al. (1997), De Haan and Ferreira (2006), Jessen and Mikosch 
(2006), or Omey and Segers (2009) for more details on regularly varying functions and max-domain of attractions. 



3 Distributional Properties of Beta-Product Convolutions 

In this section we discuss the stochastic model (jl.ip with S being beta distributed with parameters a and /?. Since 
we assume that H{0) — 0, then Hajj{Q) — 0, and H^^p possesses a positive density function ha_p given by 

KA^) = ^^^^a;"-i(J'^,p__,^,i/)(x), x£{0,Lo) (3.1) 

implying for any k GlN 

poo 

hi,k{x) = k {v-xf-^y''dH{x), xe(0,w). (3.2) 



The density function hi^k is a fc-monotone function, see Balabdaoui and Wellner (2010) for recent deep results concern 
estimation of fc-monotone functions. Conversely, any integrable fc-monotone function has representation p.2p . see 
Lemma 1 in the aforementioned paper. 

Clearly, if H possesses a density function h. then 

{Jp,P-^^p+^H){x) = {Ip,p^^^p+^h){x), X £ (0,w). 

In particular we have for some (5 G [0, 1), ?i GlV 

hi,n-s{x) ^ T{n + 1 - S){In-s,ps-„h){x), x £ {0,uj). (3.3) 

Let denote the n-fold derivative operator (we write alternatively Z^"-* instead of D^"'^ f for some diffcrentiable 
function /). If h^^^_g exist almost everywhere, utilising (|3.3p we can recover h for 5 £ (0, 1) as 

x'-hix) = ^J~\^"_^^ (/./^l"L)(x), x£iO,u:), (3.4) 



which follows by the properties of the Wcyl fractional-order integral, see Lemma [6. II and Corollary 2.1 of Fakes and 
Navarro (2007). When 6 = Lemma [HT] yields further 

h{x) = ^€(0,..) (3.5) 

which follows also from Lemma 1 in Balabdaoui and Wellner (2010). A more general result is stated in Theorem 2.1 
of the Fakes and Navarro (2007). Namely, if h exists and h £ Ii+q-^ {h £ T^-s is instead assumed therein, which is 
a misprint), then 

h{x) = (-1)" , ^^"^ _^ .T"+"-^-i(/,-i^(")(pi_„fe^.„_,-))(x), .Te(0,c.), (3.6) 
1 (a + ?i — 0) 
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provided that exists almost everywhere. When a e [0, S] formalising we arrive at: 

Theorem 3.1. Let H,Ha_^p he as above where a G [0,6], jS ~ n ^ 6,6 & [0, 1) with n G iV. Let h,haji denote the 
corresponding density functions of H and Ha^p, respectively. If h^^\ exists almost everywhere, then p.6p holds. 

Example 1. a) Consider the case a £ (0,oo) and /3 = d eJV. If exists almost everywhere and further h G Ii+q, 
then p.6p implies 

hix) = {-l)d-EM^:,-+d-iij(d)(^p^_^h^,)(^x), x€{0,u). (3.7) 
i [a + a) 

b) Suppose that a = 1/2 and (3 = d — 1/2, d <ElN . If hfj^ d~i/2 ^^^ists almost everywhere, then by Theorem 13. II 

m = (-l)'^^^'-'(/i/2^^''bi/2/ii/2,rf-i/2))(a;), xe{Q,Lo), (3.8) 
whieh reduces for /3 = 1/2 to 

h[x) = -r(l/2)(/i/2i^^i)(Pi/2/ii/2,i/2))(a;), 2;e(0,c^). (3.9) 

Example 2. Let Ha,p be the df of Ta^p^\, a Gamma random variable with positive parameters a, A and density 
function given by 

A" 

ha aix) = x"~^ exp(— Ax), x € (0, oo). 
T{a) 

Equation p.6p implies that h is the density function of T^+is^x, a Gamma random variable with parameters a + /3, A. 
If Ta+i3.\ is independent of Ba^p this means 

which is a well-known property of gamma and beta random variables, see e.g., Galambos and Simonelli (2004). 

A key fact when dealing with independent beta products is that if Bx ^ is beta distributed with positive parameters 
A = a + /3,7 being further independent of B^.p, then we have the stochastic representation 

Ba,pB\^j = Ba,fi+f. (3.10) 

The above stochastic representation is crucial for the recursive calculation of ha^fj. Since h need not always exist, it 
is of some importance to recover the df H when Ha^p is known. Utilising p.lOp this can be achieved iteratively as 
shown in our next result. 

Theorem 3.2. Let H,Hajj be two distribution functions of the random scaling model (jl.ip . If H{0) = and 
:= P > Pi > ■ ■ ■ > Pk > Pk+i := 0, fc G {0,iV} are constants such that Pi-i — /?,; G (0, l),i = 1, . . . , fc + 1, then 
there exist distribution functions Hq := H,Hi, . . . ,Hk+i = Ha^p determined iteratively by 

7J,_i(x) = r(a + /3,) ha + l3,){h^p.^.p^.iH,){x) - ( J^.,p^„_,. gz)(x)j , xe(0,w), (3.11) 

L{a + Pi-i) L ■ J 

with 6i := 1 + Pi — Pi-i. Furthermore, Hi, i = 1, . . . , fc + 1 possesses a density function hi. 

We illustrate next p. lip by two examples. 

Example 3. Consider H, Ha^p with a S (0, oo) and P = d £ IN. With Pi = d — i,i = 0, . . . ,d there exists Hi with 
density function hi,i = 1, . . . ,d such that Hq = H, Hd = Ha,p and 

H,-i{x) = ^-:E"+^-i[(a + A)a;-"-'''-^F,(x)-a--"-'5'/i,(a;) 
a + Pi I 

= H,{x)-^^^, x£{0,u;), t £{!,..., d}. (3.12) 

a + Pi 

Consequently, if h^^^^^ exists, then we can calculate h recursively by 

hi-lix) = — — h,{x) ^— 3-' ^ ^ °°)' 

a + Pi a + Pi 
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which is an alternative calculation to p.7l) . Note that if a = 1, then H can be determined explicitly by the inverse 
of the Williamson d-transform, see Proposition 3.1 in McNeil and Neslehova (2009), or Lemma 1 in Balabdaoui and 
Wellner (2010). 

Example 4. In a financial context assume that an investment R ^ H (positive) is being subjected to some deflation 
effect such that the return after a period of time (say a year) isW = RS with deflator S being uniformly distributed 
on the interval (0, 1). The fact that S = and ()3.1ip imply that H and the df H\^i of W are related by 

H{x) = Hi,i{x) - xhi^iix), yxe{0,Lu), (3.13) 
with CO € (0, go] the upper endpoint of H. Furthermore, if h^i \ exists, then almost surely in (0, uj) 

h{x) = xh^i\{x). 

Consequently, we have 



Next, if for some constant 7 



limx/ii i{x) — \iuixhi^i{x) = 0. 



xhx,\{x) ^ 11 
lim „ ' , = 7 G p, 1 
xio ffi 1 (x) 



then by Proposition 2.5 in Resnick (2007) G Further, p.l3|) implies 

lim H—h: = 1 - 7> 

x^ca Hi,i[l/X) 

which can be also written alternatively as (set R* := l/R) 

P{R*^ > x} 1 



lim 



00 P{R* > x} 1-7' 

where ^ = 1/Bi,i is a Pareto random variable with parameter 1. Karamata's Theorem (see e.g., De Haan and Ferreira 
(2006), Resnick (2007)) yields thus if 7 S [0, 1), that also H S RVy. Another proof of this fact is given in Proposition 
5.2 of Maulik and Resnick (2004). By (|3.13p a converse result can be easily established. Note in passing that since 
h{x) ~ x}i\\{x) we have h G i?Ky_i if and only if G i?Ky_2,7 



4 Lower Tail Asymptotics 

In the recent contributions Hashorva and Pakes (2010), Hashorva et al. (2010) discuss the asymptotic behaviour of 
the survival function Hafi assuming that H belongs to some max-domain of attraction of a univariate extreme value 
df. 

Hashorva and Pakes (2010) shows that R and Ha^p, belong (if so) to the same max-domain of attraction. The practical 
importance of these findings is that for insurance and finance random scaling models do not allow the deflator to 
change the max-domain of attractions of the random payment. Interesting asymptotic results for our random scaling 
models can be found in the context of Archimedean copula in Charpenticr and Sogers (2007, 2008, 2009), see Remark 
(|4.2p below. For such copulas in the aforementioned papers the asymptotics of their density generator -0 at is 
derived. The connection with our scaling model is immediate since '0 equals the survival function G^V. 

Complementing the findings of Hashorva and Pakes (2010) we focus next on the lower tail asymptotics of Ha.p, which 
boils down to determination of the min-domain of attractions of Ha p. 

When dealing with positive random variables the min-domain of attraction, for say the df H , is determined by the 
max-domain of attraction of the df of 1/i?. Since _ff is a df with lower endpoint 0, only the Frechet or the Gumbel 
max-domain of attraction for i?* is possible. The first assumption (Frechet ) is equivalent with H satisfying ()2.3|) 
with some positive index 7. Actually, this situation is simple since by Lemma 4.1 and Lemma 4.2 in Jessen and 
Mikosch (2006) if H satisfies ([Q]) with some 7 G (0,oo), and S G RVa with a G (0, 00), then 

Ha^f) G RV^^, 7* := min(7,a). (4.1) 

Note in passing that when a = 7 the asymptotic behaviour of random product follows from the well-known result of 
Embrechts and Goldie (1980). 

The beta-product convolution model is included in the above assumption since if 5 ^ ^a,^ then 

p{s<s} = (i + "(i)) r(a + nrL) ^"' ^^'^^ 
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implying 5* G RVa- 

In the next theorem we consider initially the case i?* is in the Gumbel max-domain of attraction, and then prove a 
converse asymptotic result for the Frechet case. 

Theorem 4.1. Let R ^ H be a positive random variable which is independent of S = Ba.p, and define the beta-product 
convolution df Haj3 with lower endpoint via the random scaling model We have: 

b) If 1/R has df in the Gumbel max-domain of attraction, then Ha^p G RV^. 

c) If Ha^f} satisfies (|2.3p with some 7 g (0, a), then H G RV^ and h^.p € RV^^i. 

Remark 4.2. 

1) In the setup of Archimedean copula Charpentier and Segers (2009) consider the asymptotics at of with ip the 
generator of some Archimedean copula (in the notation of McNeil and Neslehova (2009)). In the light of the findings of 
the aforementioned paper = {Hi^i)~^ , i.e., it is the inverse of the survival function of a beta-product convolution 
(a = = 1). By Proposition 2.6 (v) in Resnick (2007) S RVa, a G [0,oo] implies that i/i.i is regularly 
varying at infinity with index 1/a. For the general /c-dimensional Archimedean copula -0^^ = {Hi^k-i)~^ ,k G IV. 
Consequently, the findings of Charpentier and Sogers (2009) concern the asymptotics of the survival functions H, iJi.fe 
and that of hi,i{x) as x — > 00. Note further that the identity (|3.13p of Example 4 can also be utilised to deal with 
these functions. 

2) If a G (0, oo),/3 = d, then by (|3.7p regular variation of h holds if the same is true for D^'^\pi-aha,d)- Clearly, for 
a = 1 the latter reduces to h^a\- 

3) For any a, A G (0,oo) and S a positive scaling random variable G RVa is equivalent with S G RVaX- Conse- 
quently, our asymptotic results above apply also when is a beta random variable. 



5 Applications 

In this section we provide two applications. Motivated by the findings of Kabluchko (2010) we derive first the joint 
asymptotic independence of sample minima considering bivariate elliptical random vectors, which was shown in the 
aforementioned paper for the special case of Gaussian random vectors. 

Our second application is concerned with lower tail asymptotics for polar distributions which is also related to the 
lower tail asymptotics of aggregated risks. As mentioned in the Introduction aggregation is a central topic in various 
applications; for insurance and financial applications see e.g., Denuit ct al. (2005) and Dengcn et al. (2010). 



5.1 Asymptotics of Minima for Elliptical Samples 

When U is uniformly distributed on the unit sphere oiW& , k > 2, and A is a fc-dimensional nonsingular real matrix, 

then the random vector X = RAU is elliptically distributed. It is well-known (cf. Fang et al. (1990)) that the 
distribution function of X depends on S := AA^ but not on the matrix A itself. In view of the properties of U (cf. 
Fang et al. (1990)), if further the main diagonal of E consists of Is, i.e., E is a correlation matrix, then by Lemma 
6.1 in Berman (1983) 

X^^Xi = -Xi = RUi, l<i<k, (5.1) 
with Ui the first component of U . Furthermore, we have the stochastic representation 

A.^ ~ n i*l/2,(fc-l)/2i 

where R is independent of i3i/2.(fe-i)/2- Clearly, the random variable \Xi\ is a deflation of i? by S* = ■\/-Bi/2,(fe-i)/2- 
Next, (|4.ip and the above stochastic representation (recall also (14.2^ *) imply \Xii\ has df Q G RV^,j G (0, 1] if for 
instance H G RV^, 7 G (0, 1]. For such df Q we define constants a„, n > 1 asymptotically by 

2nP{a-^ > Xn > 0} = 1. 
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It follows that = L{l/n)n , with L G RVq a positive slowly varying function at 0, i.e., lini„_^oo L(c/n)/L{l/n) = 
1, Vc e (0, oo). For such constants we have the convergence in distribution as n — oo 

anMmA-Mi^Gj, M„i := min l<i<k, 

l<j<n 

with given by 

g^{x) = l-exp(-xT), x>0. (5.2) 

Since for u G (0,oo) small enough 

P{\Xu\<u,\Xij\<u}^0, l<i^]<k 

we have 

niO P{\Xii\<u} ' - r-l- 

Consequently, the sample minima 7Vf„ = {Mni, . . . , Mnk) has asymptotic independent components, meaning that the 
joint convergence in distribution 



holds with Ail, . ■ . , Aik independent with df G^. 

To this end, we note that the joint convergence in distribution above can be reformulated for the more general class 
of asymptotically elliptical random vectors, see Hashorva (2005) for asymptotic properties. 



5.2 Aggregation of Two Risks 

If X is a /c-dimensional elliptical random vector as above, then for given constants < i < k we have 

i<i<fe y i<j<k 

which is a well-known property for the Gaussian random vectors. Moreover, if Be{q), q € (0, 1] denotes the df of a 
Bernuli random variable assuming values —1,1, for any pair Xi , Xj ,i j we have 

(X„ X,) = (Oi, p,,Oi + y'l - PIO2) , (Oi,02) = {R,TiS,R,T2Vl - S^), S^^Bi/2,1/2, (5.3) 

with pij G (—1,1) the ijth entry of the correlation matrix S, and (Oi,02) a bivariate spherical random vector 

with positive associated random radius i?, such that i?J = R^Bi k/2-1 and Ti = Be{l/2),i = 1,2. Furthermore 
Ti, T2, i?*, S* are mutually independent. 

Since we are concerned with asymptotic results, the distribution assumption on S above can be dropped. We consider 
next the lower tail asymptotics of a bivariate polar random vector {X, Yp), p G (—1, 1) with stochastic representation 

{X, Yp) = (TiRS, pTiRS + PT2RV1 - 5"') , P ■= Vl~P^, (5.4) 

where Ti = Be{qi), qi ^ [0,1], R ^ H and S ^ G such that 

0(0) = i?(0) = 0, G(l) = l. 

As in the elliptical setup here again Ti, T2, i?, S* are assumed to be mutually independent. Since |X| = RS, then the 
lower tail asymptotics of \X\ can be established by Theorem 14.11 under asymptotic assumptions on both R and S. 
We note in passing that asymptotic properties of {X,Yp) with p E (—1,1) random are discussed in the recent 
contribution Manner and Segers (2009). Further, remark that as mentioned in Remark 14.21 we do not need to specify 
the df G, apart from the asymptotic condition G G RVa- 

We derive next the lower tail asymptotics of \Yp\ under the following additional assumption: For all positive t small 
enough 

G{p + t)-G{p-t) = Lp{t)e^, G{p + t)-G(p-t)^Lp{t)e-^, ap,ape[0,cx)), (5.5) 
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with Lp, Lp e RVo two positive functions. Condition (|5.4p can be easily checked. In the special case that G possesses 
a positive density function g continuous at p and p condition (|5.5p is satisfied with 

ap = ap = l, andLp(t) = (2 + o(l))5(p), = (2 + o(l))5(p), U 0. (5.6) 

We have now the following result. 

Proposition 5.1. Let {X,Yp)^p £ (0,1) he a bivariate random vector with stochastic representation (|5.4p . Suppose 
that 

R^HeRV^, S^GeRVa, a,7e(0,oo) 

and (|5.5|1 holds with some Oip.a-p and Lp,Lp. Assume further that when ap = ctp, then Lp(x) — cLp{x),yx > with 
some positive constant c. Then 

\X\eRV-,,, \Yp\eRV^,, (5.7) 

where 71 = min(Q!, 7) and 72 = min(7, a^, a^r). 

A simple instance for which we can apply Proposition 15 . II is when G possesses a continuous positive density function 
g. In view of (j5.6p the index 72 equals min(7, 1). The following corollary is an immediate consequence of the above 
discussion. 

Corollary 5.2. Let G, H, [X, Yp), p G (—1,1) be as in Provosition \5.1[ and let (X„, Yn), n > I be independent bivariate 
random vectors with the same df as {X,Yp). Lf G € RVi, and it possesses a continuous positive density function g, 
then we have the joint convergence in distribution 

[an min |Xj|,&„ min |Yj|) -4- {Mi,M2] n-^oo, (5.8) 

with Mi,A42 independent with c?/ ^/niin(7,i) defined in (j5.2l) and constants an,bn,n > 1 satisfying P{\X\ < l/a„} = 
■f {l^pl < ~ ^/n- for all large n. 

6 Further Results and Proofs 

Next we present first two lemmas and then proceed with the proofs of the claims in the previous sections. 
Lemma 6.1. Let /3,c be positive constants, and let h : (0, cxj) R be a given positive measurable function. 

a) If for some r E (0,oo) we have Jp x^~^\h{x)\dx < 00, then Iph is continuous at 0, i.e., 

limigh = Ink = h. (6.1) 

b) If he If} and when (3 G (0, 1) 

{y - x)^~^h{y) dy < 00, x>0 

holds for some 6 > x, then {Ifjh){x) is finite and continuous for all x > 0. 

c) If h e Ifj+c, then 

Iplch = IJph = Ifi+,h. (6.2) 

d) If the n-fold derivative D^^^h exists almost everywhere and D^"^h G Ip, then 

D'-^hph = IpD^^'^h. (6.3) 

e) For any df H with 11(0) — and upper endpoint lo G (0, 00] we have 

{Jp+i^p,_^H){x) = x{If,p^i^pH){x), x£{0,uj). (6.4) 

Proof of Lemma 16.11 By the assumption on h its Mellin transform exists, consequently utilising statement (d) 
Lemma 2.2 in Pakes and Navarro (2007) establishes the proof. Statement b) is mentioned in the Introduction of 
Hashorva and Pakes (2010). Both c) and 0?) are shown in Lemma 8.1 of the aforementioned paper. (|6.4p follows 
immediately from (|2.ip and (|2.2p . and thus the proof is complete. □ 
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Lemma 6.2. Let Ti,T2 be two random variables taking values —1,1 with P{TiT2 = —1} G (0,1] being independent 
of the scaling random variable S ^ G with G(0) = 1 — G{1) = 0. For given p G (0, 1) set Sp :~ \pTiS + /5T2 \/ 1 — S*^ | . 
If G satisfies ()5.5p . t/ien loe have 



P{Sp<u} = (l + o(l))gi,_i(H"^^M^) + (l + o(l))g-i,i(P")"^ip(w), ^iO, (6.5) 
where qij := P{Ti =i,T2= j},i,j G {-1, 1}. 

Note in passing that if G possesses a positive density function g continuous at p and p, then (j6.5p reduces to 
P{Sp<u} = il + o{l))2P{TiT2^-l}[gip)p + gip)p]u, uiO. 

Proof of Lemma [6.21 By the assumptions S G (0, 1) ahiiost surely, and Tj,j = 1,2 assumes only two values {—1, 1}. 
Hence wc may write for any u G (0, 1) small enough 



P{Sp<u} = P{Ti = l,T2 = -l}P{\pS-pVl-S^\<u} 

+P{Ti = -1,T2 = l}P{\pVl- - pS\ < u}. 
Using further the fact that S is independent of Ti , T2 we obtain 



P{\pS -pVl- S^\ <u} = P{-u< pS -py/l- <u} 

= (l + o(l)) / dG{s) 



lp-{l+o{l))pu 

il + oil)){pur^Lpiu), uiO. 



As above we have further 



P{\pVl-S^-pS\<u} = {l + o{l)){puy^Lp{u), tiiO, 
thus the result follows. □ 



Proof of Theorem 13.21 Theorem 3.3 in Hashorva and Pakes (2010) shows an iterative formula for calculating the 
survival function H when the survival function Ha. 13 is known. Our proof here is established with the same arguments 
of the aforementioned theorem utilising further p.lO|) . □ 

Proof of Theorem 14.11 a) It is well-known (see e.g., Resnick (1987)) that if R* = l/R has df in the Gumbel 
max-domain of attraction, then E{{R*y} < 00 for any (5 > 0. Hence since Ba^fs is regularly varying at with index 
a the claim follows by Breiman's Lemma (see e.g., Jessen and Mikosch (2006), Denisov and Zwart (2007), or Resnick 
(2007)). 

b) We show the proof utilising Theorem l3.2l With the notation of the aforementioned theorem there exist distribution 
functions := H, Hi, . . . , H^+i = H^ji determined iteratively by 

H,.i{x) = + A-) a+ft_, \{a + fii){Is,p-a-p,-iH,)[x) - {Js,.p_^^, Hi){x)\, V2;G(0,oo), 

with 5i := 1 + Pi — Pi-i G [0, 1) and /3o := /? > /3i > • • • > /?fc > Pk+i '■= 0. By the assumption on Ha^p for any x > 
we have (set A^.+i ;= a + Pk+i + 1) 

/>CXD 

T{5k+i){h,^,P-x,^,Hk+i){x) = / {y~xf-^--^y-'^''^-Hk+i{y)dy 

/oo 
{y - lf^^^-^y~^'^+^Hk+i{xy)/Hk+i{x) dy. 

By Karamata's Theorem (see Embrechts et al. (1997), or Resnick (2007)) 

n5u+i){h,^,P-x,^,Hk+i){x) = {l + o{l))x-"-P''Hk+i{x) {y-lY''+--^y-^>'^-+'^dy 



= (1 + o(i))..— /^-g,+,(.) ^^^^-+^);;(^^-+^ - - ^) , x;o. 

1 (Afe+1 - 7) 

Similarly we obtain 

POO 

r(<5fe+i)(J5,^,,p,_, ,ilfc+i)(x) = / [y - xY''^^-W-^''+^ dHk+i{y) 
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Jl 



(1 + o(l))7(/5,+iP_A,+ii?fe+i)(a;), X I 0. 



Consequently, 



Hk{x) 




r(a + /3fc+i + 1 - 7) 



r(a + /3fc - 7) 



i/fc+i(a;) 



(1 + 0(1)) 



r(a + /jfc+i)r(a + /3fc-7) 
r(a + /?fc)r(a + (3k+i - l) 



hence G /JK,,. Proceeding iteratively wc find tliat = H <E RV-^. 
Next, in view of (|6.6p if i7 G i?Ky, 7 G (0, 00), then as above we obtain 




dy, X 4, 



X 



estabhshing thus the claim. 



□ 



Proof of Proposition 15.11 The proof follows with the same arguments as in the proof of Theorem 14.11 applying 
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